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Lepton anomalous magnetic moments - a theory update 
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Standard Model contributions to the electron, muon, and tau lepton anomalous magnetic moments, ai — 
(gi — 2)/2, are reviewed and updated. The fine structure constant is obtained from the electron g e — 2 and used to 
refine the QED contribution to the muon g^ — 2. Recent advances in electroweak and hadronic effects on <? M — 2 
are summarized. Examples of "New Physics" probed by the a M Brookhaven experiment E821 are outlined. The 
prediction for a T is also given. 



1. Introduction and summary 

During the past few years, there has been grow- 
ing interest in leptonic anomalous magnetic mo- 
ments. On the experimental side, the extraor- 
dinary measurements of a e = {g e — 2)/2 at the 
University of Washington currently provide the 
best determination of the fine structure constant, 
a, when compared with theoretical predictions: 

a- 1 = 137.035 999 59(38)(13). (1) 

Similarly, a new effort in progress at Brookhaven 
National Lab (Experiment E821) aims to improve 
the measurement of a M by a factor of 20 or better. 
Although not competitive with a e in precision, a M 
is much more sensitive to electroweak loop effects 
as well as "New Physics" which give contributions 
~ rrif, i.e. 4xl0 4 enhancement in a M relative to 
a e . Besides being able to observe the electroweak 
loop corrections to predicted by the Standard 
Model, E821 is capable of detecting the presence 
of "New Physics." For example, supersymmetry 
loop effects can potentially provide a large con- 
tribution to a M . If a significant deviation from 
Standard Model expectations is observed, super- 
symmetry is likely to provide the leading candi- 
date explanation. 

To exploit experimental progress requires de- 
tailed calculations of the Standard Model contri- 
butions to a e and a^. QED computations at the 
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4 and 5 loop level, hadronic vacuum polarization 
and other loop effects, 1 and 2 loop electroweak 
effects, and "New Physics" contributions must be 
throughly scrutinized and refined. That effort 
challenges our theoretical tools and abilities. It 
provides an important synergism between theory 
and experiment. 



Table 1 

Standard Model predictions and experimental re- 
sults for the leptonic anomalous magnetic mo- 
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In this paper we update some of the recent 
progress in theoretical calculations of a e and a M 
and thereby provide a status report (see Table 
0). In the case of a M , we briefly describe a 
few examples of the "New Physics" sensitivity of 
E821 underway at Brookhaven. Comparison of 
the present experimental results with the Stan- 
dard Model prediction gives the following 95% 
C.L. bound on "New Physics" contributions to 
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the muon anomalous magnetic moment 

- 710 x l(T n < a™ cw Physics < 2210 x 10 -11 . (2) 

For completeness, we also give an updated pre- 
diction for a T , even though an experimental mea- 
surement of that quantity is far from current ca- 
pabilities. 



2. Electron 

To match the present experimental precision 
one has to include the following QED contribu- 
tions to a e = (g e — 2) /2 
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where B n (l, /') describe the contributions of loops 
containing lepton I' to ai, while A n contain pure 
QED contributions §-§: 
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and 
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1.1812415 
-1.5098(384) 
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5.197 x 10~ 7 
1.838 x 10" 9 
-7.3739 x 10~ 6 . 
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B 2 (e,l) describe loops with the lepton I inserted 
in the Schwinger diagram; they are calculated us- 
ing 
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B 3 (e, fi) is a sum of three groups of diagrams |1 
S 3 (e,M) = ^ (4 ' 2 »K.K)+5 (y) K/m e ) 



(7) 

They describe, respectively, diagrams with either 
two muon loops or photon corrections within a 
single muon loop; and the remaining diagrams 
with a single muon loop (they contain either a 
photon not attached to the muon loop or and 
electron loop). Finally, £?3 7 (e,/i) is a contribu- 
tion of light-by-light scattering with a muon loop 
jl2| . Their numerical values are 
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The hadronic contributions arise from vacuum 
polarization insertion in the Schwinger diagram 
frill , in the two-loop QED diagrams , and from 
the hadronic light-by-light diagram, estimated as 



-.had 



(light-by- light)m 2 /m 2 ; the result is 



had 



1.63(3) x 10 



-12 



(10) 



The electroweak contribution up to two loops is 

El 



a™ = 0.030 x 10" 12 . 



(11) 
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In total, the current Standard Model prediction 
for a e is given by 

a e = 0.5- - 0.328478444 (-) 

7T \7T/ 

(qj \ 3 / oi \ 4 

-J - 1.5098 (-J + 1.66 x 10~ 12 

(12) 

where the last term consists of the electroweak 
and hadronic contributions, whose dependence on 
a can be neglected for the purpose of a determi- 
nation at the present level of accuracy. The QED 
part of the expansion appears to converge very 
well, with alternating signs and no considerable 
growth of the coefficients. We can use the average 
of the latest experimental results for a e Jl|] , 

a^ p = 1159652188.4(4.3) x 10~ 12 , 

a e l p = 1159652187.9(4.3) x 10~ 12 , (13) 

to deduce a value of a: 

a' 1 = 137.035 999 59(38)(13) (14) 

where the first error comes from the uncertainty 
Aa° xp = 4.3 x 10~ 12 (dominated by systematic 
cavity effects), and the second from the theoreti- 
cal uncertainty in A±. 

The next most precise method of determining a 
is based on the Quantum Hall Effect (for a review 
see |l6[ ] and references therein), where one finds 


a^iqH) = 137.036 003 70(270). (15) 

We have used that value of a in obtaining the 
prediction for af 001 * in Table 0. 

Although there is currently a small discrepancy 
(1.5a) between a^ 1 (a e ) and a~ 1 (qH), the over- 
all level of agreement is truly impressive and rep- 
resents a triumph for QED. However, the good 
agreement is not a particularly severe constraint 
on "New Physics." If one assumes that such ef- 
fects contribute Aa^ cw Ph y sics ~ to 2 /A 2 , where A 
is the scale of "New Physics," then a comparison 
of a~ 1 (a e ) and a^^qH) is currently sensitive to 
A < 100 GeV. To probe the much more interest- 
ing A ~ O(TcV) region would require an order of 
magnitude improvement in a e (possibly feasible; 
see jl8| for a discussion), an improved calculation 



of A4, and a better independent direct measure- 
ment of a -1 . The last requirement could be met 
by improving a^ 1 (qH) by 2 orders of magnitude 
or perhaps more likely combining the already pre- 
cisely measured Rydberg constant with a better 
m e determination. 

3. Muon 

3.1. QED contribution 

Because of the presence of virtual electron 
loops, higher order QED contributions to a M are 
enhanced in comparison with a e . At present we 
need 5 terms of the expansion in a: 

4r-±cgy (i6) 

n=l 

with 

Ci = A x = 0.5, 

C 2 = A 2 + ax(m e /m fl ) +a 2 (m T /m fl ) 

= 0.765 857 388(44), 
C 3 = A3 + C 3 ^(e) + C 3 ^(r) + C 3 vac - poL ( e ) 

+ c vac. pol. (r)+c vac. Pol-^) 

= 24.050 509(2), 
d = Ai + 127.55(41) = 126.04(41), 
C 5 = 930(170). (17) 

where 01,2 are given in |]l9|j20|| . Taking m^/m e = 
206.768273(24) || and m T = 1777.05(26) MeV, 
m e = 0.51099907(15) MeV [| we find 

ai(m e /m p ) = 1.094258294(37), 

a 2 {m T /m ll ) = 0.000078059(23) (18) 

For the evaluation of C 2 the errors in a% j2 have 
been added in quadrature. 

In C3 we have contributions from light-by-light 
scattering diagrams with e and r loops, and vac- 
uum polarization diagrams with either e, or r, or 
both types of loops. For light-by-light we use the 
formulas of |l2] ] . For vacuum polarization we use 
p"l|], with exception of the mixed e — r diagram 
which we evaluate numerically using the kernel 
from |Q. With updated values of m^/m e and 
m T we find 

Cl J {e) = 20.9479246(7) 
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0.0021428(7) 
1.9204551(2) 
-0.0017822(4) 
0.0005276(2) 



(19) 



Even adding all errors in C3 linearly, we end up 
with the uncertainty 20 times smaller than the 
previous update |l5| . This is because of the avail- 
ability of the analytical result for A3 ||, whose 
error previously dominated. 

For C4 we use the difference between the muon 
and electron coefficients found in p2[ , and the 
latest A4 value; with errors added in quadrature 
we get 



C* 4 = 127.55(41) + At = 126.04(41). 



(20) 



The first partial evaluation of C5 was per- 
formed in p3fl , where one-loop vacuum polariza- 
tion insertions in the lowest order light-by-light 
diagram were computed, giving a contribution of 
about 570 to C5. Those are the diagrams with the 
maximal power of lnfm^/nie). However, these 
are not the only strongly enhanced diagrams. As 
had been shown in [^4j the light-by-light dia- 
grams, in which the electron loop is connected 
with 2n + 1 photons to the muon loop, contains 
7r 2 ™ lnfm^/me). The numerical coefficient of such 
term in C5 has been calculated in p5j . Some 
other diagrams, including two-loop vacuum po- 
larization insertions, were estimated in pq |; we 
adopt the value of C5 from this paper. 

Some diagrams contributing to C5 have been 
evaluated analytically in p7p^ ] . Other estimates 
of C5 have also been made, using various methods 
including renormalization group p^j29[ | and Pade 
approximation ]30fl . A class of contributions has 
even been computed to all orders f31j| . 

We now understand the ratio of the growth of 
the Ci coefficients describing the QED contribu- 
tion to the muon g — 2. In the terms calculated so 
far, the characteristic increase of the coefficients 
is of the order of 10 for each power of a /it. The 
primary reason of the growth is an extra factor of 
ln(m M /m e ) ~ 5.3 due to electron loops inserted 
in the photon propagators, further enhanced by 
combinatorial factors (this also explains why the 
higher order terms have the same sign as the first 



one; the perturbative series tries to compensate 
for the fact that we are using a too small scale, 
and therefore too small value of a). Another 
source of increase is due to the factor of ir 2 ~ 9.9, 
which accompanies each two additional photon 
rungs added to the light-by-light diagrams. This 
reasoning justifies the truncation of the perturba- 
tive series at the fifth term. Our final estimate of 
the QED contribution to is 

q qed = n6 584 705.6(2.9) x IQ- U . (21) 

The error has been estimated by linearly adding 
roughly equal contributions from the uncertainty 
in a and 64,5, and a small number estimating the 
higher order terms in the QED series for a^ ED . 

3.2. Hadronic contributions 

The bulk of hadronic contributions comes from 
vacuum polarization insertion in the Schwinger 
diagram and can be calculated using experimen- 
tal data on e + e~ — > hadrons or hadronic r decays, 
and dispersion relations (see |32| and where 
further references can be found). The most recent 
evaluation MM gives 



_ had 



(4th order) = 6924(62) x 10" 



(22) 



That study relies in part on theoretical assump- 
tions which are subject to some debate. For ex- 
ample, a preliminary update of the analysis ]33] | 
arrives at a larger error estimate of 119 x 10~ n 
II- 

In addition, one has to include the vacuum po- 
larization insertion in all two-loop QED diagrams 

us 



.had 



(6th order) = -100(6) x 10 



-11 



(23) 



and hadronic light-by-light diagram |3q| for which 



we take an average 35 1 of the values given in |37], 
|3l 



a£ ad (light-by-light) = -85(25) x 10" 
The total hadronic contribution is 



,had 



6739(67) x 10" 



(24) 



(25) 



Its error dominates the present theoretical pre- 
diction of a 

It would be very valuable to have at least an in- 
dependent estimate of the hadronic light-by-light 
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contributions from lattice QCD. Such a study 
might be undertaken in the near future ||| . Also, 
ongoing studies of e + e~ — > hadrons and hadronic 
tau decays could further reduce the theoretical 
uncertainty. A goal of ±40 x 10~ n or smaller 
is well matched to the prospectus of experiment 
E821 at Brookhaven which aims for that level of 
experimental accuracy. 

3.3. Electroweak contributions 

At the one loop level, the Standard Model pre- 
dicts 



a™(lloop) = - 



5 G^m 2 



3 8V27T 2 
\{l-As 2 w ) 2 +0 




195 x 10" 



(26) 



where G M = 1.16639(1) x 10~ 5 GeV~ 2 , M = M w 
or Mniggs , and the weak mixing angle sin 2 9\v = 
s 2 w = 1-M^-/M| = 0.224. We can safely neglect 
the O (ml/M 2 ) terms injpf). 
Two-loop corrections |45|,|15|,||,||] 



EW 



(2-loop) = -44(4) x 10" 



(27) 



decrease the electroweak contribution by about 
23%, bringing it to 



EW 



151(4) x 10 



-li 



(28) 



This decrease is mainly due to the large logarith- 
mic terms ~ ln(M§ /m 2 ) . These leading logs have 
recently been resummed to all orders in a 
They found that the leading logs in order 0(a 3 ) 
increase the EW effects by a very small amount, 
a^ w (3-loop, LL) = 0.5 x 10~ n which is within 
the error. 

The complete Standard Model prediction for 
a u is 



,SM 



- „QED 



-.had 



RW 



= 116 591596(67) x 10" 



(29) 



Currently, the combined CERN § and BNL § 
measurements give 



= 116 592 350(730) x 10" 



(30) 



It is expected that the ongoing run at BNL will 
reduce the uncertainty by a factor of three and in 
the long term ±40 x 10~ n is achievable. 

Comparing experiment and theory, one finds 



,SM 



(750 ± 733) x 10" 



(31) 



which implies good consistency; but leaves open 
the possibility of "New Physics" contributions to 
a p as large as ~ 0(1O~ 8 ). 

We subsequently consider several examples of 
"New Physics" effects on a^. To estimate the 
sensitivity of ongoing measurements and theory, 
we assume ±100 x 10~ n combined precision is 
attainable. However, if the theoretical error can 
be further reduced and the experiment proceeds 
as planned, a reduction in that uncertainty by 2-3 
appears feasible. 

3.4. "New Physics" and Electroweak Ra- 
diative Corrections 

"New Physics," beyond Standard Model ex- 
pectations, will in general give rise to addi- 
tional contributions which we collectively call 
fl Ncw Physics _ g e f ore discussing specific examples, 
we consider the electroweak radiative corrections 
to such "New Physics" effects. 

Most "New Physics" effects contribute directly 
to the dimension 5 magnetic dipole operator. In 
that case, they are subject to the same EW sup- 
pression factor as the W loop contribution to 
a^ w . From the calculation in Ref. jL5|, one finds 
a leading log suppression factor 



1 



M 



4a , 
— In 

7r m 



(32) 



where M is the characteristic "New Physics" 
scale. For M ~ 100 GeV, that factor corresponds 
to about a 6.4% reduction. 

3.5. Supersymmetry 

The supersymmetric contributions to stem 
from smuon-neutralino and sneutrino-chargino 
loops |4^-|5^1 (see Fig. 1). They can be signifi- 
cant if the supersymmetric particles are not too 
massive and if tan (3 = v-i /v\ is large. Indeed, the 
one loop effect is given in the large tan f3 limit by 



,SUSY 



57T sin 



■ tan j3 



(33) 



G 




Figure 1. Supersymmetric loop contributing to 
the muon anomalous magnetic moment. 



where m represents a typical SUSY loop mass. 
Including the EW suppression factor in Eq. (^) 
then implies 



I S , USY ~ 140 x 10- 11 



100 GeV x 



tan/3. (34) 



For large tan/3 ~ 40 scenarios, even an experi- 
mental sensitivity of ±100 x 10~ n probes in at 
the 750 GeV level and becomes competitive with 
direct high energy collider searches. 

3.6. Origin of Muon Mass 

In models where the muon mass is generated by 
quantum loops, similar loop effects will also give 
additional contributions to a M [ j53| . Under very 
general assumptions, the induced 8a „ is given by 



6a f = C JT> 



(35) 



where A is the scale of "New Physics" responsible 
for generating m M . Examples of such mechanisms 
include: extended technicolor, multi-Higgs mod- 
els, compositeness, etc. For a^ cw ph y slcs sensitiv- 
ity of ±100 x 10~ n and C ~ 0(1), A > 3 TeV is 
probed. 

3.7. Anomalous W couplings 

We generalize the jWW coupling such that the 
W boson magnetic dipole moment is given by 

e 



-(1 + K+X) 

and electric quadrupole moment by 



!w — 



2mw 



(k — A) 



(36) 



(37) 



where k = 1 and A = in the Standard Model, 
i.e. gw — 2. For non-standard couplings, one 
obtains the additional one loop contribution to 
a M given by 



a M (K, A) 



(K-l)ln 



A 2 



1 



-A 



, (38) 



where A is the high momentum cutoff required to 
give a finite result. It presumably corresponds 
to the onset of "New Physics" such as the W 
compositeness scale, or new strong dynamics. 
Electroweak effects reduce that contribution by 
roughly the suppression in Eq. ([32]). Probing a M 
at the ±100 x 10~ n level provides a sensitivity 
tok- 1| of about ±0.1 (for A - 1 TeV) (see also 
pqj). Generally, one might expect k ~ (m^/A) 2 
in theories with underlying strong dynamics at 
scale A. 

The ongoing experimental effort to improve the 
accuracy of a M measurement will improve the con- 
straints on the "New Physics" scenarios we men- 
tioned above, as well as on other theories and 
phenomena, such as a general two-Higgs doublet 
model |59tl , leptoquarks ] |60| , or four fermion con- 
tact interactions [[n) (see also Ip^l ). 

4. Tau 

The anomalous magnetic moment of the tau, 
a T , is predicted to be 



QED 



had 



EW 



where 

a QED 



,EW 



1.1732 x 10~ 3 , 
3.2(4) x 10~ 6 , 
4.7 x 10 -7 . 



(39) 

(40) 
(41) 
(42) 



We have updated those contributions to incorpo- 
rate m r = 1.77705(26) GeV, a theoretical esti- 
mate of a^ ad from perturbative QCD 0] which 
is averaged with the result given quoted in 163], 
and two loop EW effects [|l5| which suppress af™ 
by about 15%. In total, one finds 



a T = 1.1769(4) x 10" 



(43) 



Experiments are currently not sensitive enough to 
measure a T . They can, however, indirectly bound 
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an anomalously large a T due to "New Physics." 
For example, e + e~ — » t + t~ cross-section mea- 
surements imply Ba] 



|a T | < 0.02. 



(44) 



Also, at this meeting L. Taylor B reported on 
analysis of e + e~ — > t + t~7 by L3 and OPAL 
which gave a similar constraint (see Table p. 
If "New Physics" effects in a T are of the form 
m^/A 2 , then one would need to extend such con- 
straints to \a T \ < 0.0003 to probe the interesting 
regime A > 100 GeV. 
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